Non-monotonous short-time decay of the Loschmidt echo 
in quasi-one-dimensional systems 
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We study the short-time stability of quantum dynamics in quasi-one-dimensional systems with 
respect to small localized perturbations of the potential. To this end, we address, analytically and 
numerically, the decay of the Loschmidt echo (LE) during times short compared to the Ehrenfest 
time. We find that the LE is generally a non-monotonous function of time and exhibits strongly 
pronounced minima and maxima at the instants of time when the corresponding classical particle 
traverses the perturbation region. We also show that, under general conditions, the envelope decay 
of the LE is well approximated by a Gaussian, and we derive explicit analytical formulas for the 
corresponding decay time. Finally, we demonstrate that the observed non-monotonicity of the LE 
decay is only pertinent to one-dimensional (and, more generally, quasi-one-dimensional systems), 
and that the short-time decay of the LE can be monotonous in higher number of dimensions. 

PACS numbers: 03.65.-w, 03.65.Sq, 05.45.Mt 



I. INTRODUCTION 

Nearly a quarter of a century has passed since Peres 
[l[ made an important advance in understanding the ori- 
gin of irreversibility in quantum theory. He pioneered a 
framework in which dynamical instabilities were studied 
in terms of small external perturbations of the Hamil- 
ton operator without resorting to the classical concepts 
of mixing and coarse graining. In particular, Peres con- 
vincingly demonstrated that a change in the time evolu- 
tion of a quantum system caused by the perturbation is 
largely determined by whether the system exhibits regu- 
lar or chaotic behavior in the classical limit. 

The quantity addressed by Peres, currently known as 
the Loschmidt echo (LE) in the quantum chaos commu- 
nity and the fidelity in the field of quantum information, 
is defined as [l| 



M(i) = (*o|exp(iff 2 ^) ex PH#i i A)l*o/ 



(1) 



It is the squared overlap of the initial state l^o) and 
the state obtained by propagating |\&o) though time t 
under the original, unperturbed Hamiltonian H±, and 
then through time —t under the perturbed Hamilto- 
nian H2- In other words, the LE, M(t), is a mea- 
sure of the time-reversibility of a system's dynamics un- 
der "imperfect conditions" . In an alternative interpre- 
tation, the LE characterizes the distance between two 
states in the Hilbert space, I'Fi) = exp(— iHit/H)\^o) 
and 1^2) = exp(— i^t/fij^o), both obtained from the 
same initial state, |\&o)j m the course of the quantum 
evolution through time t under two different Hamiltoni- 
ans, Hi and H2 respectively. By construction, M(0) = 1 
assuming the initial state is normalized. Typically, M{t) 
decreases (or "decays" ) with increasing time, and the pre- 
cise functional form of the decay is determined by the 
nature of the Hamiltonians Hi and H2, as well as by the 
initial state. 



Interest to the subject of the LE revived a decade ago. 
This was largely due to the discovery of a perturbation 
independent regime [1] of the LE decay in systems that 
exhibit chaotic dynamics in the classical limit. In this 
regime, known as the Lyapunov regime, the LE decays 
exponentially in time, M(t) ~ exp(— At), with the decay 
rate A equal to the average Lyapunov exponent of the 
corresponding classical system. The Lyapunov regime 
can be considered as a clear example of how classical 
chaos manifests itself in a quantum mechanical world. 

Roughly speaking, for chaotic systems under the action 
of global Hamiltonian perturbations, i.e., perturbations 
affecting a dominant part of the system's phase space, 
one identifies three distinct decay regimes of the LE: the 
Gaussian perturbative regime [H, 0| for "weak" perturba- 
tions, the exponential Fermi-golden-rule (FGR) regime 
for "moderate" perturbation strengths, and the ex- 
ponential Lyapunov regime 0, @ for "strong" perturba- 
tions. The above three regimes serve as a basis for the 
theory of the LE decay for global perturbation. The full 
theory however is much more subtle and the interested 
reader is referred to two comprehensive review articles, 
Refs. 0,1- 

The LE due to global perturbations has long been a 
subject of numerous experimental studies. Most notably 
these are experiments in nuclear magnetic resonance |9|- 
ll~2l ]. quantum optics cold atoms fl3 - [l8j . supercon- 
ductivit y 1131. microwave cavities |20l - [23j . and elastody- 
namics [24|, [25[ . 

In recent years, Hamiltonian perturbations of a new 
kind, namely perturbations that are local in phase space, 
have been addressed in the context of the LE decay both 
theoretically and experimentally [26l - [30| . In a semiclassi- 
cal picture, a local perturbation is concentrated in a small 
region of the system's phase space, so that the length of 
a typical trajectory between any two successive visits of 
the perturbation region is much larger than the system's 
size. In strongly chaotic systems the phase space extent 



2 



of a local perturbation can be characterized by an "es- 
cape" rate r~* denned as the rate at which trajectories 
of the corresponding classical system visit the pertur- 
bation region. The semiclassical analysis of the LE for 
times longer than the Ehrenfest time, i.e., longer than 
the time it takes for an initially localized wave packet 
to explore the available phase space, has revealed an ex- 
ponential decay regime, M(t) ~ exp(— nt), in which the 
decay rate k is a non-monotonous function of the per- 
turbation strength (2(| [Hj]. In particular, for sufficiently 
weak perturbations the FGR regime is recovered in which 
k grows quadratically with the perturbations strength. 
However, in the limit of strong perturbations k saturates 
at a perturbation independent value 2t~^ corresponding 
to the so-called escape-rate regime. The crossover from 
the FGR to the escape-rate regime is non-monotonous, 
and k exhibits well pronounced oscillations as a function 
of the perturbations strength. These oscillations have 
been then confirmed in numerical experiments with per- 
turbed cat maps |29j and, more recently, in laboratory 
experiments with microwave cavities [301 ] . Finally, the 
limit of point-like perturbations, for which r osc — > oo, 
was addressed in Ref. [27|- There, the LE was shown to 
decay algebraically with time, M(t) ~ t~ 2 . 

The existing theory of the LE decay from local per- 
turbations is only applicable to times long compared to 
the Ehrenfest time. However, in certain cases the short 
time decay of the echo, during which the unperturbed 
and perturbed quantum states can be described by local- 
ized wave packets, might be of significant importance for 
understanding results of laboratory experiments. For in- 
stance, in echo spectroscopy experiments with cold atoms 
trapped inside an optical billiard [l5|-[l7j one typically 
observes the echo decay for times as short as only few 
free flight times of the corresponding classical billiard. 
Such time scales can be short compared to the Ehrenfest 
time depending on a parameter choice. It is the objec- 
tive of this paper to address the short-time decay of the 
LE due to local Hamiltonian perturbations. More specif- 
ically, we make the first step in this direction by per- 
forming a detailed study of the LE in one-dimensional 
(and, more generally, quasi-one-dimensional) systems in 
the presence of localized perturbations of the system's po- 
tential. We show that in such systems M(t) is typically 
a non-monotonous function exhibiting well pronounced 
minima and maxima at times when the particle traverses 
the perturbation region. We also show that in closed 
systems the envelope of M(t) can be well approximated 
by a Gaussian, exp [— (t/r) 2 ] , with the decay time r ex- 
plicitly expressible in terms of system parameters. All 
results presented in this paper concern the LE decay in 
"clean" systems and imply no averaging over initial states 
or Hamiltonian perturbations. 

The analysis presented in this paper only concerns the 
short-time decay of the LE in conservative quasi-one- 
dimensional, and therefore essentially integrable, sys- 
tems. The case of the LE decay in classically chaotic 
systems for times shorter than the Ehrenfest time (and 



in the presence of global perturbations) was addressed 
in Ref. [3l|. There, for a "typical" localized initial state 
the LE was shown to exhibit the double-exponential 
initial decay, M(t) ~ exp(— constant x e 2A ') with A being 
the Lyapunov exponent. 

The paper is organized as follows. In Section |H] we 
provide the analysis of the LE decay in one-dimensional 
systems based on full numerical solution of the time- 
dependent Schrodinger equation and on the thawed 
Gaussian approximation (TGA) in its standard and mod- 
ified, average potential formulations. Details on the av- 
erage potential TGA are deferred to Appendix [AJ Par- 
ticular examples treated in Section |TT] include a free par- 
ticle (Sec. Ill Ap . particle on a ring (Sec. IIIBI) . harmonic 
( Sec. Ill C~Tj) and anharmonic (Sec. Ill C~2|) oscillators. Sec- 
tion IIIII demonstrates disappearance of non-monotonous 
features of the LE decay as a quasi-one-dimensional sys- 
tem is transformed into a substantially two-dimensional 
system. In Section [IV] we give a discussion of our results 
and make concluding remarks. 



II. LOSCHMIDT ECHO IN ONE DIMENSION 



A. Free particle 



■7-H 




FIG. 1: (Color online) One-dimensional quantum particle in 
the presence of a localized potential barrier. See the text for 
discussion. 



The essential physics of the short-time decay of the 
LE for local potential perturbations can already be re- 
vealed by considering the simplest physical system - a 
free particle. In this scenario, the unperturbed Hamil- 
tonian, Hi = p 2 /2m, describes one-dimensional motion 
of a free particle of mass m; in the coordinate repre- 
sentation, the momentum operator p = —ihd/dq. The 
perturbed Hamiltonian, H2 = H\ + V(q), represents a 
particle interacting with a potential "barrier" V(q) lo- 
calized to an interval on the position axis, see Fig.[TJ For 
concreteness, we assume V{q) to have a single maximum 
(or minimum) at q — a. (As it will become clear from 
the following presentation, our results can be easily gen- 
eralized to perturbation potentials of arbitrary shape.) 
We also label the characteristic extent of the potential 
by 7, see Fig. [1] Finally, we set the initial state to be a 
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Gaussian wave packet 
i 

1 



*o(g) 



exp ( ^p (q - qo) ^ — ) ( 2 ) 



that represents a quantum particle with the average po- 
sition coordinate q and momentum p Q . The initial wave 
packet dispersion is quantified by a. We further assume 
that the initial state, \E , o(?)j nas zer0 overlap with the 
perturbation potential, i.e., 



7,cr < I, 



(3) 



where I = \a — qo\ is the characteristic distance between 
the particle and the perturbation region, see Fig. [T] We 
also assume that the initial state is well localized in the 
momentum space, so that 



p a > h, 



(4) 



and that the potential V(q) constitutes a classically small 
perturbation: 



\V(q)\<ZE, 



(5) 



where E — p\j2m is the total energy of the corresponding 
classical particle. Together Eqs. Q and (JSJ imply that 
the probability density gets almost perfectly transmitted 
over the barrier in the course of the time evolution, and 
the effect of quantum reflections can be neglected. 

Most of semiclassical studies of the LE decay rely on 
the so-called "dephasing representation" (DR) [H, [H| , 
in which the LE amplitude (\E f 2 (^) l^i (^)) is expressed as 
an interference integral over trajectories of the unper- 
turbed system, modulated by a "dephasing" factor due 
to the perturbation, with initial phase space coordinates 
weighed by the Wigner function of the state |^o). Our 
analysis however is based on a different analytical ap- 
proach, namely on the "thawed Gaussian" approximation 
(TGA) |34l - [3a ] in its standard and extended versions, see 
Appendix [A] and the discussion below. The TGA, unlike 
the DR, involves only a single trajectory transporting the 
wave packet center through phase space in the course of 
its time evolution. This allows for greater analytical flex- 
ibility at the expense of a reduction in accuracy. For most 
of our purposes however the TGA proves to be sufficiently 
reliable both qualitatively and quantitatively. 

According to the TGA both the unperturbed and per- 
turbed wave functions preserve their Gaussian form in 
the course of their time evolution: 



(ma. 



x exp -a 3 {q - qj) 



1j) 



(6) 



where j = 1,2 labels the unperturbed and perturbed 
wave functions respectively, and 5i stands for the real 
part. Here, qj(t), Pj(t), fa(t) are real-valued and a.j(t) 



complex- valued functions of time that par ametrize the 
wave functions. In the "standard" TGA |34l - l3q ] the time 
evolution of these parameters is governed by 



Pi = -W;(qj), 



2ih 



2h 



<(<&)> 



(7) 
(8) 

(9) 



and^ =Pj/2rn—Wj(qj)—(h 2 /m)Uaj. Here Wj(q) is the 
corresponding (unperturbed for j — I and perturbed for 
j = 2) potential, and the prime denotes differentiation 
with respect to q. The system of four real, first order 
ordinary differential equations given by Eqs. (J7HH]) is then 
solved with the initial conditions qj — qo, pj — po, and 
aj 1 = 2(T 2 at t = 0. In the "average potential" TGA, 
see Appendix El Eqs. (J7HSJ) are replaced by 



qj =pj/m, 

Pj = -W l j [q j -^cx j ), 



(10) 
(11) 
(12) 



with the average potential Wj defined as 

Wj(q;s) = (2s/tt)5 [ dxW j (x)e- 2s ^- q ^ . (13) 



As before the prime denotes partial differentiation with 
respect to q. 

In the free particle case the unperturbed wave function 
^i(q;t) given by the TGA, Eq. flB}, is exact and the 
parameters qi, pi, and a\ evolve in time according to 



9i 
Pi 



qo - 
Po ■ 



p t/ m , 



2(a 2 + iht/m) . 



(14) 
(15) 
(16) 



and 4>\ — (po/2m)t— (h/2) arctan(fe/mc7 2 ). 

Time evolution of the parameters (72, P2, «2 character- 
izing the perturbed wave function ^(g; t) are determined 
from a set of three ordinary differential equations (JTHH]) 
in the "standard TGA" and by equations (fTUl - [T3"|) in the 
"average potential TGA" with j = 2 and W 2 (q) = V(q). 
Here we note that the LE defined by Eq. (pj is unaffected 
by the global phases (f>i and fa in the regime that al- 
lows one to approximate the unperturbed and perturbed 
states by the simple Gaussian wave packets given by 
Eq. ([B]). The phase <$> of an individual Gaussian wave 
packet would only be physically relevant if the initial 
state ^0 was a superposition of two or more Gaussian 
wave packets, or if one was interested in the LE amplitude 

2{t)\^ i{t)) rather than in M(t). However, considera- 
tions of these kinds go beyond the scope of the current 
study. 
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Substituting the two time-dependent Gaussian wave 
packets, given by Eq. (J^J) with j = 1 and 2, into Eq. ([T]) 
we obtain for the LE 



M (t) =-r— — n- ex P 



\a>x + a 2 
x [gft(aia 2 (ai + a 2 )*)Ag 2 



+ 3(aia 2 ) 



AqAp 
ft 



(ai + a 2 J 



Ap 2 
' 4ft 2 



, (17) 



where Ag = q 2 — <?i and Ap = p 2 — pi, asterisk de- 
notes complex conjugation, and 3 stands for the imagi- 
nary part. Equation (|17jl. along with equations describ- 
ing time evolution of the parameters qj , pj , and ay with 
j = 1,2, provides the main framework for our analytical 
study of the LE. 

We also perform a numerical analysis of the problem 
of the LE decay. To this end we adopt a method of ex- 
pending the quantum propagator exp(— ifft/ft) in terms 
of Chebyshev polynomials of the Hamiltonian H . For a 
detailed discussion of the method see Ref. [37j and ref- 
erences within. Hereinafter we refer to results of the nu- 
merical solution of the LE decay problem as to "exact" 
results as opposed to approximate ones obtained by using 
the standard and average potential TGA. 
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FIG. 2: (Color online) Top figure: The LE, M(t), for a 
free particle with the perturbation potential V(q) given by 
Eq. QB) . The horizontal dashed line represents the LE satu- 
ration value M(oo) given by Eq. (|28p . Bottom figure: Pertur- 
bation potential as seen by the unperturbed classical particle, 
i.e. V(qi(t)). The system is characterized by po = 50, a — 0.1, 
a — qo = 1, 7 = 0.3, and Vb = 150. 

Figure [5] shows the LE decay for the case of a pertur- 
bation potential given by the Gaussian barrier 



V(q) = V exp 



(q-af 



(18) 



The upper part of the figure displays M(t). Here, the 
thick blue curve shows a result of the full, quantum- 
mechanical computation of the LE defined by Eq. ((T|), 
i.e., the exact LE decay. The green curve is obtained 
by using the standard TGA, i.e., Eq. (fT7f with the time 
evolution of parameters q 2 , p 2 , and a 2 determined from 
Eqs. (J7HS1). The red curve represents the LE decay ob- 
tained in the average potential TGA by using Eqs. (|10f - 
1131) for the time evolution of <7 2 , p 2 , and a 2 . The initial 
wave packet, Eq. ©, is specified by q — — 1, p = 50, 
and <t = 0.1 [43j]. The perturbation potential is given by 
Eq. (US) with a = 0, 7 = 0.3, and Vb = 150. The lower 
part of Fig. [2] shows the function V{q\{t)) 1 which has the 
meaning of the perturbation potential as "seen" by the 
unperturbed classical particle. 

Figure [2] prompts an interesting observation: the LE 
develops a minimum at the time the unperturbed clas- 
sical particle passes the extremum of the perturbation 
potential, and a maximum at the time the particle ex- 
its the perturbation region. In fact, such behavior is 
generic. This can be shown by imposing the approxi- 
mation a 2 (£) ~ o:\(t), which simplifies Eq. (|17[) to 



M(t) = exp 



1 ( 2 2 AqAp Ap 2 

'iiad |ai1 Aq +5ai T + 4ft 2 " 



Then, substituting Eq. (fT6|) into Eq. (fl"9| we obtain 



(19) 



M(t) = exp 



2a 



Aq Ap) -- 1 



(20) 

The function M(t), given by Eq. pO)) . possesses two ex- 
trema at time instants determined by 



3^ = ° 



Aq = 



ht 



(minimum) , (21) 



— Ap (maximum) . (22) 
m 



We now utilize the standard TGA, Eqs. © and ® 
with j = 2 and W 2 (q) — V(q), to express the posi- 
tion and momentum shifts, Aq and Ap, in terms of the 
perturbation V(q). To the leading order in Vo/E, with 
Vb = max„ |V(g)|, we obtain 



1 



Aq(t) = -- I dt'V( qi (t')) , 
Pa Jo 

TTi 

Ap (t) = V( qi (t)). 

Po 



(23) 
(24) 



Substitution of Eq. ((24|) into Eq. jSTJ shows that M(t) 
exhibits a minimum at time t a = m(a — qo)/po = ml/\po\, 
with a being the extremum point of the perturbation 
potential, see Fig. Q] It is the time instant at which the 
unperturbed classical particle passes the extremum point 
of the perturbation potential, i.e., q\{t a ) = a. Similarly, 
the LE exhibits a maximum at time % = m(b — qo)/po 
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that, in accordance with Eqs. (|22l) . (|23l) and (|24|) . satisfies 



B. Particle on a ring 



U dt'V( gi (t')) = (t b + (m °"^ )2 ) ^(gi(t 6 )) . (25) 

In the limit t 3> ma 2 /h, or equivalently / 3> po<r 2 /h, 
Eq. (E5J) simplifies to 



dg^(g) = |6-gb|V(6). 



(26) 



Equation f|26[) has a simple geometrical interpretation as 
depicted in Fig.[TJ the area under the curve V(q) from qo 
to b equals the area of a rectangle with the base | b — qo | 
and height V(b). It is clear from this construction that 
the LE exhibits a maximum when the classical unper- 
turbed particle leaves the perturbation region. 

The time instant t = t&, at which M(t) exhibits a max- 
imum, has the following physical significance. According 
to Eq. (|22l) . and in the limit I 3> poa 2 /h, it satisfies 
Aq(tf,) — &.p(tb)tb/m. This means that if the perturbed 
wave packet, ^2(9;^), is propagated backward through 
time — tf, under the unperturbed (free-particle) Hamilto- 
nian Hi then the resulting wave packet is centered at go- 
In other words, the states |* ) and | e ^i*f>M e -^2 VM^) 
have the same expectation values of the position opera- 
tor (while generally different expectation values of the 
momentum operator). 

Finally, at long times, t 3> ml/po, one has Aq = 
— {m/po) dqV(q) and Ap = 0. Thus, in accordance 
with Eq. (|2"U|l. the LE saturates at the value 



M(oo) = cxp 



1 

aE 



dqV(q) 



(27) 



For the case of a perturbation potential given by Eq. ([T8")) 
the saturation value is 



M(oo) = exp [-(7r/8)( 7 /a) 2 (K)/^) 2 



(28) 



In Fig. [2] the LE saturation value, M(oo), calculated in 
accordance with Eq. ([28} is shown by a horizontal dashed 
line. We note that the quantitative agreement between 
the LE saturation value obtained by numerically solv- 
ing the Schrodinger equation and the saturation value 
given by Eqs. (j2"7| and (|28p improves as the perturbation 
strength Vo is reduced. This is consistent with the fact 
that Eqs. (|27|) and (|28|) are valid only in the leading order 
of V /E. 

As it will become clear from the following sections the 
observed non-monotonicity of M(t) is a generic feature 
of the LE decay in one-dimensional (and, more generally, 
quasi-one-dimensional) systems. We now turn our at- 
tention to closed systems and address the envelope of 
M(t). The analysis presented here is valid for short 
times, for which the wave function stays localized and 
can be well approximated by a Gaussian wave packet. In 
other words, we restrict our discussion to t <C £e with 
being the Ehrenfest time. 



As the fist example of a closed one-dimensional system 
we consider a quantum particle moving on a ring. In 
this set-up the wave function of the unperturbed system, 
^>i(q;t), is defined on an interval < q < L, and sat- 
isfies the time-dependent Schrodinger equation with the 
Hamiltonian Hi = p 2 /2m and periodic boundary condi- 
tion ^i(L;t) = >3>i(0;t). The Hamiltonian of the per- 
turbed system is H 2 = H 1 + V(q), where V{q) is defined 
on the same g-interval with V(L) = V(0). As before, we 
assume V(q) to be localized in a small region of size 7 
around a single extremum point a, see Fig. [TJ The wave 
function of the perturbed system, "^(g; t), satisfies the 
same periodic boundary conditions, ^2(L;t) — ^ 2 (0;t). 
Finally, as in Section III A[ we assume the validity of con- 
ditions given by Eqs. ([5HS]). 

As stated above, we only consider the wave packet dy- 
namics for times shorter than the Ehrenfest time, t <§; £e- 
The latter is the time that it takes for the wave packet 
dispersion, initially equal to a, to become comparable to 
the size of the system, L. (We note, that Eq. ([3]) implies 
a <C L.) According to the uncertainty principle the ve- 
locity uncertainty of the particle is Av ~ h/ma. Thus, 
for the case of the free motion on a ring, the Ehrenfest 
time can be estimated as £e ~ L/ Av ~ maL/h. 




0.025 



FIG. 3: (Color online) Top figure: The LE, M(t), for a par- 
ticle on a ring with the perturbation potential V(q) given 
by Eq. (|18|) . The dashed line represents the envelope decay, 
M(t), calculated in accordance with Eqs. (|29(l and (|31[). Bot- 
tom figure: Perturbation potential as seen by the unperturbed 
particle, i.e. V(qi(t)). The system is characterized by L = 1, 
po = 200, a = 0.1, a - q = 0.3, 7 = 0.1, and V = 1500. 

Since the dispersion of the wave packet and the ex- 
tent of the perturbation potential are small compared 
to the ring circumference L, the problem in question is 
effectively equivalent to that of a free unbounded mo- 
tion in the presence of a periodic perturbation potential 
J2n=-oo V(q + nL). Therefore, according to Section \U~K[ 
we expect the LE curve to exhibit a sequence of minima 
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at times i m i n ,fc = m(a—qo+kL) /po, with k S No, at which 
the unperturbed classical particle traveling a constant ve- 
locity po /m passes the extrema of the periodic perturba- 
tion potential. (For concreteness but without loss of gen- 
erality, we consider the case < qo < a < L and po > 0.) 
Indeed, such a non-monotonous time decay of the LE is 
found in perfect agreement with the results of the fully 
numerical (exact) and semianalytical (TGA) solutions of 
the problem, see Fig- [3j We also note here that, just as in 
Fig. [21 M(t) in Fig. [3] exhibits a sequence of local maxima 
at times i+ ax fc ~ m{a - q + 7 + kL)/p and t m ^ k+1 ~ 

?7i(a — qo — 7 + (fc + l)i) /po, with k G No. These time in- 
stants are solutions of Eq. (|25|) for tj, with V(q) replaced 
by the periodic potential J2ti=-oo V{l + n L)- 

In order to find the envelope function M(t) of the LE 
decay curve we consider M at times t n = nT, where n £ 
N, and T = mL/po is the period of the unperturbed mo- 
tion. According to Eqs. (|20"|) and (l23|) we have M(t n ) = 



exp(~Aq 2 /2a 2 ) 
leads to 



with Aq 



-(t n / Po L) dqV(q). 



M(t)=cxp[-(t/T) 2 ] 
with t = r r i nc , and 



Tring = V2\p \aL 



dqV(q) 



This 



(29) 



(30) 



For the case of a Gaussian perturbation potential given 
by Eq. ([T5]) the expression for the decay time reduces to 



T 



aE 
ttj\V \ ' 



(31) 



where E = p^/2m is the total energy of the system. 

The dashed curve in Fig. [3] is the envelope function 
M(t) calculated from Eqs. (|2"9")l and (j3"Tj) for the system 
specified by L = 1, po = 200, a = 0.1, a — qo = 0.3, 
7 = 0.1, and V = 1500. It is evident that M(t) ac- 
curately describes the envelope of the exact LE decay 
represented by the thick blue curve. It is also interest- 
ing to compare the quality of the approximate solutions. 
The standard TGA (green curve) reasonably well approx- 
imates M(t) around its maxima, while completely failing 
to capture the function close to its minima. On the other 
hand, the average potential TGA (red curve) well de- 
scribes both maxima and minima. This approximation 
however turns out to be only limited to short times, up 
to 3 full cycles of the particle around the ring for the par- 
ticular set of parameters, after which the approximation 
becomes unstable. 



C. Particle in a well 

We now consider a quantum particle trapped inside a 
one-dimensional potential well, so that Hi = p 2 /2m + 



E 



U(q) 



Po 



q- qo 



a 



FIG. 4: (Color online) A quantum particle in a one- 
dimensional potential well. The unperturbed potential is 
given by U(q), and the perturbed potential by U(q) + V(q). 
The total energy E of the system is such that both unper- 
turbed and perturbed classical motions have the same turning 
points q- and q+. 



U{q) and H2 = p 2 /2m + U(q) + V(q), where U(q) speci- 
fies the potential well, and V(q) represents a perturbation 
potential localized around q = a. As before, we assume 
validity of the conditions given by Eqs. ([3] [5]) with the 
total energy of the classical particle E — p 2 , /2m + U(qo). 
We also assume E > U(a) + V(a), so that both unper- 
turbed and perturbed classical motions have the same 
turning points g_ and q + , see Fig. [4] 

As we will see below the LE decay for the system under 
consideration is non-monotonous. Deep local mimima 
surrounded by (generally less pronounced) local maxima 
occur when the classical unperturbed particle traverses 
the perturbation region. The physics underlying this 
non-monotonicity is essentially the same as in the cases of 
a free particle and a particle on a ring, see Sections III Al 
and III B\ so we directly proceed to the analysis of the 
envelope function M(t) of the LE decay M(t). 

We begin by considering the period T + AT of the 
classical oscillatory motion of the perturbed system: 



T + AT = V2m 




dq 



^E - U(q) - V(q) 



- J dqV{q){E-U{q)Y\ (32) 



where 



/,_ ^E - U(q) 



(33) 



is the period the unperturbed classical motion. Equa- 
tion ((32)) holds to the second order in V/ (E — U). Then, 
assuming that (E — U(q)) is approximately constant in 
a 7-interval about q — a, we obtain 



AT=^(E-U(a)Y 



+ OO 



dqV(q) . 



(34) 



Here we note that if the perturbation potential satis- 
fies J dqV(q) = then the expansion in Eq. (|32|) has 
to be terminated at the next order resulting in AT ~ 
JdqV 2 (q). 
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The envelope function M(t) of the LE decay curve can 
be calculated using the approximation given by Eq. (|19|) . 
To this end we consider the distance between the un- 
perturbed and perturbed trajectories in phase space, 
A<7 = <72 — <7i and Ap = p 2 — pi , at times t n — nT with 
n £ N. Since T is the period of the unperturbed systems 
we have qi(t n ) = qo and pi(t n ) = po. The trajectory of 
the perturbed system however is "delayed" by the time 
nAT with respect to the unperturbed trajectory. This, 
for small AT, can be written as q2(t n ) = qi(t n — nAT) 
and P2(t n ) — Pi(tn — nAT). Expanding these equations 
to the leading order in AT/T we obtain 



. . po AT 

Aq(t n ) = — t n , 

m 1 

AT 

Ap(t n ) = U'(q )—t n , 



(35) 
(36) 



where the prime denotes the derivative. Then, a substi- 
tution of Eqs. $55$ and into Eq. (JTSJ) yields 



M(t n ) = exp 



-X(tr. 



AT 



(r 



-t, 



(37) 



with 

X = 



foot! 



(Kl^) 
V mJ 



PiU'{qi) , fU'( qi ) 



mh 



2h 



(38) 

keeping in mind that q\ = qi(t), p\ — pi(t), u\ = a\{t), 
and qi(t n ) = q , Pi(t n ) = pa- 

It is now interesting to observe that x is a constant of 
the motion defined by Eqs. ([?]-[§]) with j — 1 and W\(q) = 
U(q). Indeed, it is straightforward to verify that 



^=0. 
dt 



(39) 



Then, a replacement of x(i n ) in Eq. (|37p by its value at 
time t = 0, 



X 



1 (Po \ 2 1 ( vU'jqo) 



2 \ma 



(40) 



yields the Gaussian LE envelope M(t) given by Eq. ([29 
with the decay time r = t wc \\ and 



T i T 



/ Po 



<rU'(q ) 



(41) 

Below we consider some particular potentials U(q) and 
compare the analytical prediction given by Eq. (|4ip with 
the numerical, full quantum-mechanical solution of the 
LE decay problem. 



1. Harmonic oscillator 

As our first example we considering motion of a particle 
in a harmonic potential well for which U(q) — mw 2 g 2 /2. 



In this case the time evolution of the unperturbed system 
is governed by [36| 



qi = qo cos lot - 



Po 



sinwt . 



Pi — Po cos ut — mioqo sin tot , 



/mw\ h cos tot + imuja 2 sin cot 



\ 2h J ih sin cot + mtoa 2 cos cot 



(42) 
(43) 

(44) 



Equations and (|4"2"H4"4")) give the exact quantum dy- 
namics of the system with the Hamiltonian H\ and initial 
state of Eq. ([2]). As before, the time dependent param- 
eters 92, P2, and «2 characterizing the perturbed wave 
packet are calculated using the standard TGA, Eqs. (JTJ- 
IHl), or the average potential TGA, Eqs. ([TUIITBl) . with 
W 2 (q) = U{q) + V(q). Then, the TGA approximation 
of the LE is computed from Eq. (fT7|) . and the envelope 
function, M(t), from Eqs. ([29]) and (gl]) with T = 2-k/lo 
and AT given by Eq. 
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FIG. 5: (Color online) Top figure: the LE, M(i), for a par- 
ticle in a parabolic well with the perturbation potential V(q) 
given by Eq. (|18[) . The dashed line represents the envelope 
decay, M(t), calculated in accordance with Eqs. (|29l) and (I41|) . 
Bottom figure: the perturbation potential as seen by the un- 
perturbed particle, i.e. V(qi(t)). The system is characterized 
by lo = 60, go = —1, po = 10, a = 0.1, a = 0, 7 = 0.1, 
and Vq = 200. Note that the LE decay obtained by the aver- 
age potential TGA (red curve) is essentially indistinguishable 
from the exact result (blue thick curve) for most of the time 
range. 

Figure [5] shows the time decay of the LE for a quan- 
tum particle in the harmonic potential well. The un- 
perturbed system is defined by u = 60, qo = — 1, 
Po = 10, a = 0.1. The perturbation potential is given 
by Eq. HB]) with a = 0, 7 = 0.1, and V = 200. As 
before, the thick blue curve represents M(t) obtained 
by numerically solving the time dependent Schrddinger 
equation for the perturbed system. (The exact time evo- 
lution of the unperturbed system is given by Eqs. © 
and (|42| - f4"4"l) .') The green curve corresponds to the stan- 
dard TGA for the perturbed system, i.e., Eqs. (J7HS]) with 
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W2 (g) = muj 2 q 2 /2 + V(q). It is clear from the figure that 
the standard TGA fails to reproduce the minima of the 
exact LE curve. On the other hand, the average poten- 
tial TGA result, computed according to Eqs. (|10HT5|) and 
represented by the red curve, provides an excellent ap- 
proximation of the exact result. The dashed line shows 
the envelope decay, M(t), predicted by Eqs. and 
(l4"Tj) . A reasonable agreement between the theoretical 
and numerical results is apparent. 



and W2(q) = U(q) + V(q). It well approximates the ex- 
act LE decay curve for up to two full oscillations of the 
unperturbed system. The dashed line shows the envelope 
dec ay, M(t ), predicted by Eqs. and (|H|). As in Sec- 
tion III C 11 one observes a reasonable agreement between 
the theoretical and numerical results. 



III. BEYOND ONE DIMENSION 



2. Anharmonic oscillator 

We now consider, as our unperturbed system, a quan- 
tum particle moving in an anharmonic potential well of 
the form U{q) = mui 2 q 2 /2 + eq 3 . Here, the computa- 
tional procedure of the LE decay is essentially the same 
as in Section III C II with the only difference that the os- 
cillation period T is obtained by numerically evaluating 
the integral in Eq. 
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FIG. 6: (Color online) Top figure: the LE, M(t), for a particle 
in an anharmonic well with the perturbation potential V(q) 
given by Eq. (I18[) . The dashed line represents the envelope 
decay, M(t), calculated in accordance with Eqs. (|29[) and (|41|l . 
Bottom figure: the perturbation potential as seen by the un- 
perturbed particle, i.e. V(qi(t)). The system is characterized 
by cj = 60, e = -400, q = -1, po = 10, a = 0.1, a = 0, 
7 = 0.1, and Vo = 200. 

Figure [6] shows the time decay of the LE for a quantum 
particle in the anharmonic potential well. The unper- 
turbed system is defined by U) — 60, e = —400, qo = —1, 
Po = 10, a = 0.1. The perturbation potential is given 
by Eq. (I3J) with a = 0, 7 = 0.1, and V = 200. The 
thick blue curve represents the exact LE decay. The 
green curve corresponds to the standard TGA for the 
unperturbed and perturbed systems, i.e., Eqs. ([TH^J) with 
W x {q) = U{q) and W 2 (q) = U(q) + V(q). Once again, the 
standard TGA fails to reproduce the minima of the exact 
LE curve. The red curve shows the result of the average 
potential TGA, i.e., Eqs. I|IDHI5|) with Wi(q) = U(q) 



In this section we argue that the reported non- 
monotonicity of the short-time decay of the LE is per- 
tinent mainly to one-dimensional (or, more generally, 
quasi-one-dimensional) systems. The disappearance of 
the non-monotonicity of the LE decay in systems with 
two or higher number of dimensions can be understood 
as follows. In one-dimensional systems an initially Gaus- 
sian wave packet traverses the perturbation region and 
approximately preserves its Gaussian shape provided the 
strength of the perturbation potential is small compared 
to the energy of the particle. In this case the TGA pro- 
vides a natural basis for the LE analysis and the theory 
developed in Section |TT] holds. In higher number of di- 
mensions however it is often energetically preferable for 
the wave packet to split in parts and circumvent the per- 
turbation region. The wave packet deforms to avoid the 
perturbation region and can no longer be approximated 
by a Gaussian. The overlap between the unperturbed 
(approximately Gaussian) and perturbed (substantially 
non-Gaussian) wave packets decays rapidly and mono- 
tonically during the time that the system interacts with 
the perturbation. 

To illustrate this point we consider the time decay of 
the LE for a free quantum particle in two dimensions 
under the action of a Gaussian perturbation potential. 
The initial state of the particle is given by 



*o(a:,y) 



■ exp ( ^p (x - x ) 



(x - x ) 2 + (y- y f 
2a 2 



(45) 



Here, (xo,yo) is the center of the wave packet, po the 
magnitude of the momentum pointing along the £-axis, 
and a its dispersion. The perturbation potential is de- 
fined as 



V(x, y) = Vq exp - 



(x - a) 2 



r 



T_ 

8 2 



(46) 



with 7 and S quantifying the extent of the perturbation 
region in the x- and y-directions respectively. As in Sec- 
tion |TT] we assume the validity of the conditions given by 
Eqs. (J3H5])- 

Figure [7] shows the time decay of the LE, M(t), 
obtained by numerically solving the time-dependent 
Schrodinger equation with the initial state given 
by Eq. (|4"5l) and the Hamiltonian operators Hi = 
(-h 2 /2m)(d 2 /dx 2 + d 2 /dy 2 ) and H 2 = Hi + V(x,y), 
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t 

FIG. 7: (Color online) The LE decay for a free particle in 
two dimensions under the action of the perturbation potential 
given by Eq. ((46}. The system is characterized by po = 50, 
a = 0.1, a - qo = 1, Vo = 150, and 7 = 0.3. The M(t) 
curves correspond to seven different values of the y-extent of 
the perturbation region: from bottom to top 5 = 0.1, 0.2, 0.3, 
0.35, 0.4, 0.5, and 00. 



where V(x,y) is defined by Eq. (|46p . The system pa- 
rameters are po = 50, a = 0.1, a — qo = 1, Vo = 150, 
and 7 = 0.3. Six different values of the y-axis extent of 
the perturbation region are considered, S — 0.1, 0.2, 0.3, 
0.35, 0.4, 0.5, as well as the case 6 — > oo corresponding 
to the one-dimensional problem analyzed in Section [II Al 
It is clear from the figure that in quasi-one-dimensional 
systems, i.e., when a <C S so that the y-component of 
the classical force exerted on the particle by the pertur- 
bation potential is negligible, M(t) exhibits minima and 
maxima as the particle traverses the perturbation region. 
On the other hand, in substantially two-dimensional sys- 
tems, i.e., for 5 < a, a fast monotonous decay of the LE 
is observed. 

In order to support our qualitative explanation for the 
disappearance of the non-monotonicity in the LE decay 
for small values of 5 we depict in Fig.[5]the time evolution 
of the unperturbed and perturbed wave functions for the 
system characterized by po — 50, a = 0.1, a — qo = 1, 
Vo = 150, 7 = 0.3, 6 = 0.1. This set of parameters corre- 
sponds to the LE decay represented by the red (most bot- 
tom) curve in Fig. Here, a snapshot at the top shows 
the probability distribution of the initial state. The left 
column shows the probability distribution of the unper- 
turbed wave packet at regular time intervals. The cor- 
responding snapshots in the right column give the prob- 
ability distribution of the perturbed wave function. A 
white ellipse in the right column shows the characteristic 
extent of the perturbation potential: the ellipse is cen- 
tered at the maximum of the Gaussian potential, and its 
major and minor semi-axes equal to 7 = 0.3 and 5 = 0.1 
respectively. The figure demonstrates how the probabil- 
ity distribution of the perturbed wave function changes 




FIG. 8: (Color online) The time evolution of the unperturbed 
(left column) and perturbed (right column) wave function. 
The top snapshot shows the initial state. The white ellipse in 
the snapshots of the right column is centered at the position 
of the perturbation maximum, and has its major and minor 
semi-axes equal to 7 = 0.3 and 8 = 0.1 respectively. The 
system is characterized by po = 50, a = 0.1, a — qo = 1, Vo = 
150, 7 = 0.3, S — 0.1, and exhibits a LE decay represented by 
the red (most bottom) curve in Fig. [7] 

from a simple Gaussian to a double-peak function as the 
quantum particle circumvents the perturbation region. 
Clearly, this process can not be adequately described by 
the TGA. It is also natural to expect the LE overlap to 
decay monotonically and faster than in the limit of large 
S, in which the perturbed wave function does not sig- 
nificantly change its functional form and only gets dis- 
placed in phase space with respect to the unperturbed 
wave packet. 

The qualitative change of the LE decay shown in Fig. [7] 
can be linked to an effective decrease of the Ehrenfest 
time is of the perturbed systems as the perturbation size 
d decreases: the observed splitting of the wave packet 
indicates that <e ~ m(a — qo)/po for small values of 8. 

It is interesting to note that the phenomenon of wave 
packet splitting is observed already for relatively weak 
perturbation potentials. For example, the right column 
in Fig. [5] corresponds to the energy of the classical par- 
ticle, Pq/2> being more than 8 times higher than the 
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strength of the potential, V$. At this point, it is not 
clear whether the splitting can be adequately explained 
in terms of deflection of the corresponding classical tra- 
jectories or whether quantum interference effects play a 
crucial role. It can be speculated that the wave packet 
splitting is intimately related to the phenomenon of co- 
herent electron flow branching in two-dimensional quan- 
tum wells with weak disordered background potentials 
[38M40I ] . However, a proper investigation is yet to be per- 
formed. 

As a final remark we make the following counter- 
intuitive observation. By increasing the parameter S, 
while keeping Vq fixed, one increases the overall pres- 
ence of the perturbation, in the sense that the integral 
J dxdy V{x, y) becomes larger. However, at the same 
time, the larger 8 the higher the LE saturation value, see 
Fig [JJ and therefore the smaller the effect of the pertur- 
bation on the time evolution of the wave function. The 
reason for that is that it is energetically preferable for 
a wave function to avoid small perturbation regions by 
changing its shape. This results in small values of the LE 
overlap. On the other hand, spatially extended perturba- 
tion regions are traversed by the wave function without 
any significant shape change resulting in large values of 
the LE overlap. 



on a ring, Eq. (|30)) . and a particle inside a smooth poten- 
tial well of an arbitrary shape, Eq. Plj) . The analytical 
formulas have been given convincing numerical support. 

All results presented in this paper pertain to the short- 
time decay of the LE in "clean" quantum systems and 
assume no averaging over initial states or Hamiltonian 
perturbations. In this respect, the findings may prove 
valuable in echo experiments which imply no or minimal 
averaging intrinsically. Thus, for example, in echo ex- 
periments with ultra-cold atoms one could try to use an 
experimentally observed decay function M(t) to recon- 
struct the localized perturbation potential. 
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IV. CONCLUSIONS 

We have studied, analytically and numerically, the de- 
cay of the Loschmidt echo (LE) from local Hamiltonian 
perturbations for times short compared to the Ehren- 
fest time. During those times the quantum states of 
the unperturbed and perturbed systems can be efficiently 
approximated by localized Gaussian wave packets. The 
latter evolve in accordance with equations of the thawed 
Gaussian approximation (TGA). We have analyzed these 
equations and deduced several important properties of 
the LE decay in quasi-one-dimensional systems. 

Firstly, we have shown that the LE is generally a non- 
monotonous function of time. More specifically, M(t) 
exhibits strongly pronounced minima and maxima at the 
instants of time when the corresponding classical particle 
traverses the perturbation region. The minima of M(t) 
are reached when the particle passes the peak (or bot- 
tom) of the perturbation potential, whereas the maxima 
generally occur at the times the particle enters or exits 
the perturbation region. We have also demonstrated that 
while the observed non-monotonicity of the LE decay is 
generic in one-dimensional systems the short-time decay 
can be monotonous in systems with higher number of 
dimensions. 

Secondly, using the TGA we have analyzed the enve- 
lope decay of the LE and shown it to be well approxi- 
mated by a Gaussian, M(t) — exp [— (t/r) 2 ] . The decay 
time t is expressible in terms of parameters of the initial 
state, system's Hamiltonian, and perturbation potential. 
We have given explicit formulas for the cases of a particle 



Appendix A: Thawed Gaussian approximations 

The celebrated thawed Gaussian approximation 
(TGA), originally introduced by Heller [3J], is a semi- 
classical technique for propagation of localized Gaussian 
wave packets in time without having to solve the full 
Schrodinger equation. The approximation is based on 
the idea that, at least for short times, the expectation 
values of the position and momentum operators, qj and 
Pj respectively, evolve according to the classical, Hamil- 
ton equations of motion. Thus, in the TGA the evolving 
wave packet is "guided" by a single real phase-space tra- 
jectory described by the Hamilton equations. The TGA 
is one of the simplest approximations to the full Van 
Vleck-Gutzwiller semiclassical propagator |4lj . 

Below we begin with a brief discussion of the standard 
TGA following Ref. [H| • A comprehensive review of the 
subject can be found in Ref. [35|]. We then introduce 
an extended version of the approximation - the average 
potential TGA - that proves more reliable for potentials 
with rapid spatial variations. 

The central assumption of the TGA is that an ini- 
tially Gaussian wave packet preserves its Gaussian form, 
i$>j(q;t) as given by Eq. ©, throughout the time evolu- 
tion under a Hamiltonian Hj — p 2 /2m+Wj(q). The wave 
packet is parametrized by one complex- valued - otj(t) - 
and three real- valued - Qj(t), Pj(t), and <j)j(t) - func- 
tions, which are determined by substituting tyj into the 
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Schrodinger equation: 



iha 



2ft_ 
m 



i — ocj ) \<i <h) + Pj (q - qj) 



2ihaj [ qj 
2m 



Pj 



W 3 (q) 



qj) + 4>j - 
ih 5id, 



4 fta,- 



PjQj 

* i (?;t) = 0. (Al) 



Here dots represent differentiation with respect to time. 
Clearly, Eq. (|A1|) can not be satisfied for a general poten- 
tial Wj (q). However, an approximate solution of Eq. (|A1[) 
can be found by expanding the potential Wj{q) into 
the power series about q = qj to the second order in 
(q — qj), and then separately comparing terms of the 
same order. This procedure yields Eqs. (JTHSJ) together 
with fa = p]/2m - Wjiqj) - (h 2 /m)^ aj . 

The standard TGA, Eqs. (THHI), relies on the quadratic 
approximation of the potential Wj (q) about q = qj, and 
is therefore limited to potentials varying slowly on the 
scale given by the spatial extent of the wave packet. In 
order to lift this limitation we modify the TGA method 
as discussed below. 

The problem of finding the "best" functions otj(t), 
qj(t), Pj(t), and (f>j(t) approximating the time evolution 
of the wave packet allows for an alternative approach. 
If *j (q; t) , given by Eq. ^ , were a true solution of the 
time-dependent Schrodinger equation with the Hamilto- 
nian Hj = p 2 /2m + Wj(q) then the expectation value of 
any (generally time-dependent) operator Oj would sat- 
isfy 



■x(*i\O j \9 j ) = -(* j \[H j ,O j ]\* 3 , 



dt 



(A2) 

with [-,-] denoting the commutator. For a general po- 
tential Wj, however, ^j(q;t) given by Eq. ^ is not a 



true solution of the Schrodinger equation, and, therefore, 
Eq. (IA2|) can only be satisfied by four linearly indepen- 
dent Hermitian operators Oj. A choice of these four op- 
erators uniquely defines the functions aj(t), qj(t), and 
Pj(t). (Equation (|A2|) is unaffected by the global phase 
4>j ■ This phase however is of no importance to the prob- 
lem of the Loschmidt echo decay addressed in this paper.) 
Although there is no unique choice of the four indepen- 
dent Hermitian operators Oj we proceed with the seem- 
ingly most natural ones: q, p, q 2 , andp 2 . Straightforward 
calculations yield Eqs. (fT0HT3]) for the time evolution of 
qj, pj, and ctj. Obviously, the requirement that ^j(q;t) 



satisfies Eq. (|A2[) for Oj = q, p, q , p guaranties the 
resulting approximation, which we here term the aver- 
age potential TGA, to predict the phase space dynamics 
of the wave packet's center and dispersion in the best 
possible way. 



Here two final remarks are in order. Firstly, since 
W^qy^aj) = (^W^g)^) Eqs. (10]) and (HJ) are 
nothing but the statement of the Ehrenfest theorem (3(| . 
Secondly, it can be readily shown that in the case of the 
potential Wj{q) being a second order polynomial in q 
Eqs. dlHlJ) are identical with Eqs. (fT0HT2|) . In this case 
*5>j(q;t) constitutes the true solution of the Schrodinger 
equation. 

The idea of exploiting average potentials for propa- 
gation of Gaussian wave packets has been previously 
discussed in the chemical physics literature, e.g., see 
Ref. [42| . These approaches however deal with poten- 
tials averaged with respect to Gaussian wave functions 
with "frozen" exponents (such as ^j(q;t) in Eq. ((HJ) with 
ctj(t) = ckq) and are generally different from the average 
potential TGA presented above. It would be interest- 
ing to compare these different methods in terms of their 
precision and computational efficiency. 
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